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ABSTRACT 


This  report  describes  experimental  and  theoretical  investigations  of  the 

particle  emission  produced  in  a  laser-surface  interaction.  Measurements 

have  been  made  on  the  angular  distribution  of  ions  emitted  from  a  tungsten 

target  irradiated  by  a  ruby  laser  at  flux  densities  around  10  to  20  megawatts/ 

2 

cm  .  In  the  plane  perpendicular  to  the  plane  containing  the  incident  laser 

radiation  and  normal  to  the  target  surface  the  number  of  ions  emitted  follows 
2 

a  cos  9  distribution  where  9  is  the  angle  from  the  normal  to  the  target.  In 

the  plane  containing  the  laser  radiation  and  the  target  normal,  most 

ions  appear  to  be  emitted  within  approximately  15c  of  the  normal  to  the 

target.  The  energies  of  the  ions  are  also  maximum  in  the  direction  normal 

to  the  target  and  fall  off  as  the  angle  from  the  normal  increases.  Calculations 

on  a  model  for  heating  an  alkali  ion  plasma  by  absorption  of  laser  radiation 

g 

indicate  that  runaway  heating  can  occur  at  flux  densities  around  10  watts/ 

2  i  8  3 

cm  with  particle  densities  around  10  /cm  .  The  magnitudes  and  time 

histories  of  shock  waves  produced  in  nickel  and  carbon  targets  by  absorption 

of  laser  radiation  have  also  been  calculated. 
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SECTION  I 


INTRODUCTION 

The  investigations  described  in  this  report  are  continuations  of  work  which 

have  been  carried  out  previously.  These  earlier  investigations  have  been 

concerned  with  particle  emission  produced  in  laser-surface  interactions. 

(1  2  3  4  5) 

This  work  has  been  described  in  detail  in  previous  reports.  '  ’  ’  ' 

Some  knowledge  of  the  background  of  this  work,  as  discussed  in  the  previous 

reports,  is  necessary  for  an  understanding  of  the  methods  and  aims  of  the 

investigations  described  in  the  present  report.  A  time-of-flight  spectrometer 

has  been  used  to  measure  ion  emission  produced  by  absorption  of  ruby  laser 

2 

radiation  with  power  densities  of  the  order  of  50  MW/ cm  .  The  ions  were 
found  to  be  mainly  alkali  metals  with  energies  of  the  order  of  several  hun¬ 
dred  electron  volts.  The  neutral  molecule  emission  produced  under  similar 
conditions  has  been  studied  with  a  quadrupole  mass  spectrometer  and  found 
to  consist  of  thermally  desorbed  gases  such  as  Hg,  CO,  and  C02  along  with 
pulses  of  high  energy  neutral  molecules  having  energies  of  the  order  of 
100  eV. 

In  this  report  we  describe  our  extension  of  these  measurements.  The  bulk 
of  the  experimental  work  reported  here  has  been  carried  out  in  an  inter¬ 
action  chamber  with  collectors  around  the  periphery  of  the  hemisphere 
which  has  the  laser-illuminated  target  at  its  center.  Measurements  of  the 
angular  distribution  indicate  that  there  is  a  maximum  in  both  the  number 
of  ions  and  the  ion  energy  in  a  direction  normal  to  the  target  surface,  and 
that  both  these  quantities  fall  off  as  the  angle  from  the  normal  increases. 
Moreover,  if  one  looks  in  the  plane  defined  by  the  direction  of  the  laser 
beam  and  the  normal  to  the  target  surface,  one  obtains  a  different  angular 
distribution  than  for  the  plane  perpendicular  to  this. 
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As  we  shall  see  later,  there  have  been  considerable  problems  with 
noise  in  obtaining  this  data.  The  essential  criterion  for  obtaining  good 
data  on  this  system  appears  to  be  the  use  of  relatively  low  laser 
flux  densities.  This  means  that  this  system  may  not  be  so  useful 
for  obtaining  data  at  high  laser  flux  densities  as  we  had  at  one  time 
hoped. 

On  the  theoretical  side  we  have  carried  out  an  analysis  of  the  model 

in  which  the  plasma  was  heated  by  absorption  of  laser  radiation  in  a 

plasma  of  alkali  ions  with  little  coupling  of  energy  to  the  desorbed 

neutral  ga.ses.  This  picture  contained  many  physically  appealing 

features  and  was  in  accord  with  our  experimental  data.  However,  the 

computer  analysis  indicates  that  heating  takes  place  only  at  laser  flux 

densities  and  at  alkali  ion  densities  which  are  uncomfortably  high. 

Therefore,  the  results  of  this  analysis  are  essentially  negative  and  it 

appears  that  invoking  this  mechanism  is  not  sufficient  to  explain 

heating  of  the  blowoff  material  at  laser  flux  densities  below  100  megawatts/ 
2 

cm  . 

We  have  also  extended  calculation  of  the  pressure  pulses  produced  by 
thermomechamcal  shock  waves  induced  ip  the  target  material  by 
absorption  of  the  laser  radiation.  Using  the  PUFF  computer  code  we 
have  calculated  the  time  history  of  the  stress  profiles  for  nickel  and 
carbon  targets. 

Finally,  we  have  calculated  steady- state  heating  for  thin  metallic 
targets.  This  topic  has  become  of  interest  with  the  advent  of  high 
power  carbon  dioxide  lasers.  The  results  indicate  that  for  a  proper 
choice  of  target  material  the  C02  laser  can  be  a  valuable  tool  in 
studying  the  laser  surface  interaction. 
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EXPERIMENTAL  WORK 

The  bulk  of  the  experimental  work  in  this  report  period  has  been  performed 
in  the  hemispherical  interaction  chamber  which  has  collectors  arranged 
around  the  periphery  of  the  hermisphere  and  a  target  at  the  center.  The 
only  target  material  used  to  date  has  been  tungsten.  A  considerable  amount 
of  data  has  been  collected  using  this  system.  Measurements  have  been 
obtained  for  the  angular  distribution  of  the  ions  emitted  from  the  tungsten 
target  irradiated  by  high  power  pulses  of  Q- switched  ruby  laser  radiation. 

A.  NOISE  IN  SYSTEM 

There  have  been  a  considerable  number  of  problems  in  the  operation  of  the 
system.  The  main  problem  has  been  that  of  noise.  The  data  obtained  from 
the  collectors  has  frequently  been  noisy.  A  typical  example  is  shown  in 
Figure  II- 1.  It  shows  two  photographs  of  high  frequency  oscillations  near 
the  beginning  of  pulses  from  the  output  of  the  collector  located  along  the 
normal  to  the  target  surface  and  from  the  hemisphere.  These  oscillations 
often  mask  the  form  of  the  output  signal  and  blot  out  portions  of  it.  The 
persistent  presence  of  the  high  frequency  oscillations  has  been  very 
troublesome  in  obtaining  good  data  using  this  system. 

The  noise  problem  is  particularly  apparent  when  the  laser  double  pulses 
This  is  shown  in  Figure  II - 2  which  shows  two  photographs  of  the  output 
when  the  laser  emitted  a  second  pulse.  In  the  upper  photograph  the 
second  pulse  occurred  about  30  microseconds  after  the  first  pulse,  and  in 
the  lower  photograph  about  20  microseconds  later.  The  emission  from  the 
first  pulse  is  relatively  noise  free;  then  at  the  beginning  of  the  emission 
produced  by  the  second  laser  pulse  there  are  high  frequency  oscillations 
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Figure  II-l  -  Noise  in  Signals  from  Interaction  Chamber 

A.  2  Us/ cm.  Top  trace  0.  02V/ cm  from  collector  along 
normal  to  target.  Bottom  trace  lOV/cm  from  hemis¬ 
phere. 

B.  0.02V/ cm,  5  |is/ cm.  Top  trace  from  collector  at  1  5° 
from  normal  to  target.  Bottom  trace  from  collector 
along  normal  to  target. 
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Figure  II- 2  -  Noise  in  Signals  from  Interaction  Chamber  Due  to  Laser 

Double  Pulse 

A.  5  Us/ cm.  Top  trace  0,  005V/ cm  from  collector  at  45° 
from  normal  to  target.  Bottom  trace  0.  02V/ cm  from 
collector  along  normal  to  target. 

B.  5us/cm.  Top  trace  0.  005V/ cm  from  collector  at  1  5° 
from  normal  to  target.  Bottom  trace  0.  02V/cm  from 
collector  along  normal  to  target. 
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Since  tne  second  laser  pulse  was  smaller  than  the  first  and  produced  less 
particle  emission,  it  is  possible  that  the  noise  originates  in  interaction 
of  tne  radiation  with  the  blowoff  material  so  that  the  system  is  particularly 
prone  to  noise  when  there  has  been  particle  emission  produced  a  short  time 
before  a  laser  pulse  comes  in. 

We  have  tried  a  number  of  techniques  for  noise  suppression.  One  of  these 
is  shown  in  Figure  11-3  which  shows  a  schematic  diagram  of  the  hemisphericai 
interaction  chamber  with  the  target  and  the  collector.  Ions  traveling  from 
the  target  to  the  hemisphere  pass  through  the  grid  to  the  collector  giving 
rise  to  a  voltage  signal  at  point  A  on  the  oscilloscope  preamplifier  If 
noise  results  from  a  common  mode  or  ground  loop  or  in  pickup  in  the 
oscilloscope  leads,  it  could  be  eliminated  by  a  differential  measuring 
technique.  We  put  the  inputs  A  and  -B  from  the  oscilloscope  preamp  at 
point  1  and  adjust  the  gain  of  A  and  -B  for  zero  output  from  the  sum  of  A 
and  •  B  using  short  leads  coiled  around  each  other.  We  then  move  -  B  to 
point  2  and  take  the  measurements.  W;e  find  that  tins  technique  makes  no 
difference  in  the  output  results  and  conclude  that  the  noise  is  not  a  spurious 
effect  caused  by  pickup  in  the  electronic  circuitry.  It  may  originate  in 
oscillations  induced  by  movement  of  charge  in  the  target  collector 
hemisphere  system.  At  present  we  do  not  understand  well  the  origin  ol 
this  noise. 

The  noise  can  be  eliminated  by  operation  of  the  laser  at  a  relatively  low 

power  so  that  the  laser  flux  density  at  the  position  of  the  target  is  less  than 

2 

approximately  20  megawatts/ cm  .  At  this  level  we  obtain  clean  output 
pulses  from  the  co.'lector  such  as  are  shown  in  early  portions  of  the 
oscilloscope  traces  m  Figure  11  2  Whenever  the  laser  flux  density  is 
kept  below  this  level  good  data  is  obtained  from  the  collectors  except 
in  the  case  whei  e  the  laser  double  pulses  and  the  emission  produced  by 
the  second  pulse  is  noisy.  However,  even  when  the  laser  double  pulses 
data  can  be  obtained  from  the  early  portion  output  of  the  collectoi  since 
the  interesting  phenomena  occur  before  the  second  laser  pulse. 
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Figure  II - 3  -  Schematic  Diagram  for  Noise  Suppression 


Since  the  oscillations  are  of  much  higher  frequency  than  the  frequencies  of 
the  desired  signals,  filtering  techniques  may  improve  the  situation. 

This  limitation  on  operating  flux  density  means  that  the  usefulness  of  the 
mteraction  chamber  may  be  limited  and  that  we  could  use  it  for  only 
certain  ranges  of  laser  power.  Above  that  power,  high  frequency  high 
amplitude  oscillations  which  block  out  a  portion  of  the  signal  make  the 
obtaining  of  good  data  difficult.  Thus,  we  may  not  be  able  to  use  this  inter¬ 
action  chamber  to  extend  the  operating  range  of  laser  powers  in  this  study 
as  we  had  hoped  to  do. 

B.  POTENTIALS  IN  SYSTEM 

The  effect  of  applying  a  potential  to  the  target  was  investigated.  We  held 
the  target  at  a  positive  potential  with  respect  to  the  hemisphere.  1.  e.  ,  an 
ion  accelerating  field.  Some  typical  results  are  shown  in  Figure  11-4. 

The  top  photograph  was  taken  at  0  accelerating  potential,  the  middle  pnoto 
graph  of  270  volts,  and  the  bottom  photograph  at  540  volts.  In  each  of  these 
photographs  the  upper  trace  is  from  the  collector  along  the  normal  to  the 
target  and  the  bottom  trace  is  from  the  hemisphere.  Both  were  taken 
across  50  ohm  resistors  and  fed  into  the  oscilloscope  The  sweep  speeds 
of  all  these  photographs  are  five  microseconds  per  cm.  On  the  top  photo 
graph  the  vertical  sensitivity  of  the  top  trace  is  0.  1  volts/ cm,  the  lower- 
trace  .  1  volt/ cm.  In  the  middle  and  bottom  photographs  the  top  traces  are 
0.  1  volts/cm,  and  the  bottom  traces  are  10  volts/ cm.  The  top  photograph 
shows  the  behavior  that  one  expects  from  the  system.  In  the  lower  trace 
one  sees  displacement  current  at  the  hemisphere  beginning  essenfally  at 
the  time  of  the  laser  pulse  and  rising  to  a  peak  at  the  time  the  charge  begins 
to  arrive  at  the  position  of  the  hemisphere  as  shown  by  the  break  in  the  tr  ace 
from  the  collector.  The  relative  shapes  of  these  two  traces  are  easy  to 
understand.  However,  the  relative  amplitudes  are  somewhat  puzzling  m 
that  the  solid  angle  subtended  by  the  collector  is  much  smaller  (  /-v/  ]  0  ’  ) 
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Figure  II- 4  -  Effect  of  Potential  Applied  to  Target.  Sweep  speed: 

5  us/ cm.  Upper  traces  from  collector  along  normal 
to  target.  Bottom  traces  from  hemisphere.  All  upper 
traces  are  0.  I  V/ cm.  In  the  top  photograph  the  bottom 
trace  is  0.  IV/ cm;  in  the  middle  and  bottom  photographs 
the  lower  traces  are  lOV/cm. 


than  the  solid  angle  subtended  by  the  hemisphere,  but  the  amplitudes  of  the 
signals  are  different  bv  only  a  factor  of  3.  Part  of  this  is  no  doubt  due  to 
the  emission  being  peaked  in  the  direction  normal  to  the  surface  however 
this  effect  cannot  explain  all  of  the  relative  differences. 

When  270  volts  are  applied  to  the  target,  some  rather  surprising  results  are 
obtained.  The  results,  as  shown  in  the  middle  photograph,  are  reproducible 
and  can  be  obtained  when  one  goes  back  to  the  same  condition  day  after  day. 
The  signal  arriving  at  the  hemisphere  increased  by  two  orders  of  magnitude 
and  the  shape  is  changed.  The  little  bump  at  the  beginning  may  possibly  be 
a  remnant  of  the  original  shape  as  shown  in  the  top  photograph.  However 
the  main  current  arriving  at  the  hemisphere  appears  to  be  delayed  by 
approximately  the  drift  time  of  the  plasma  to  the  hemisphere.  It  lias  a 
shape  very  Similar  to  that  shown  for  the  pulse  shape  at  the  collector  in 
the  top  photograph.  Moreover,  the  pulse  shape  at  the  collector  has  become 
rather  odd.  It  shows  a  rapid  decrease  after  it  reaches  the  peak.  Some 
irregular  bumps  follow  the  peak. 

When  the  voltage  is  raised  to  540  volts,  the  pulse  shape  at  the  collector 
again  becomes  regular  and  has  approximately  the  same  shape  as  when  no 
voltage  is  applied  to  the  target.  Increasing  the  target  voltage  to  540  volts 
does  not  appreciably  change  the  shape  of  the  pulse  at  the  hemisphere  as 
compared  to  its  shape  at  2  70  volts.  This  also  is  rather  surprising  since 
one  would  expect  the  charge  to  pass  across  the  gap  more  quickly  as  the 
accelerating  voltage  is  raised. 

One  might  be  able  to  explain  the  behavior  of  the  pulse  shape  at  the  hemisphere 
by  assuming  a  plasma  that  drifts  most  of  the  way  across  the  gap  between 
target  and  hemisphere  without  charge  separation.  At  higher  voltages  it  is 
pulled  apart  after  it  has  expanded  to  a  certain  point  so  that  the  displacement 
current  effects  do  not  show  up  in  the  bottom  two  photographs  of  Figure  1.1  4. 
This  also  would  explain  the  large  increase  in  the  pulse  size  at  the  hemisphei  e 
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as  one  applies  an  accelerating  voltage  for  ions,  and  it  would  explain  the  small 
dependence  of  the  time  of  the  arrival  of  this  pulse  on  the  accelerating  voltage. 
However,  if  this  were  the  case,  one  would  expect  the  pulse  size  at  the  collector 
to  be  increased  considerably  also,  whereas  the  collector  pulse  size  is  almost 
independent  of  the  voltage  applied  to  the  target.  Moreover,  there  is  no  satis¬ 
factory  explanation  regarding  the  anomalous  behavior  of  the  collector  trace 
at  270  volts.  At  this  stage  there  are  a  range  of  phenomena  in  the  interaction 
chamber  which  we  cannot  adequately  interpret. 

We  investigated  the  field  profiles  in  the  chamber  as  voltage  was  applied  to 
the  target.  Equipotent ial  surfaces  in  the  angular  distribution  chamber  were 
determined  using  an  analog  method.  The  electrode  configuration  of  the 
chamber  was  drawn  on  a  sheet  of  conducting  paper  using  silver  paste, 
and  appropriate  voltages  were  applied  to  the  different  electrodes.  One  of 
the  two  terminals  of  a  high  impedance  probe,  such  as  an  oscilloscope,  was 
tied  to  a  point  while  the  other  one  was  moved  in  the  region  between  the 
electrodes.  The  equipotential  lines  are  then  obtained  by  following  the  lines 
along  which  nulls  in  the  potential  difference  between  the  two  probes  are 
located.  Such  equipotential  surfaces  are  of  use  in  interpreting  the  angular 
distribution  measurements  that  are  made  with  a  potential  difference  between 
the  target  and  other  electrodes  in  the  chamber. 

The  equipotential  surfaces  were  determined  for  two  different  potential 
configurations.  The  first  configuration  is  shown  in  Figure  II- 5.  In  this 
case  the  target  is  held  at  a  positive  potential  of  +10  volts,  while  the 
hemisphere  and  outer  shell  of  the  chamber  are  held  at  ground  potential. 

It  is  obvious  from  the  figure  that,  except  in  the  immediate  neighborhood 
of  the  target,  the  electric  field  is  radial  at  all  points  of  interest  in  the 
chamber.  The  second  configuration  is  one  in  which  the  target  and  outer 
shell  of  the  chamber  are  held  at  ground  potential  while  the  hemisphere  is 
held  at  a  potential  of  -10  volts.  The  resulting  equipotential  surfaces  are 
shown  in  Figure  11-6.  In  this  case  it  is  obvious  that  there  are  large  tangential 
components  of  electric  field  even  at  points  far  removed  from  the  target. 
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Potential  Distribution  in  Interaction  Chamber  With  Target  at  riO  Volt 
and  all  Other  Electrodes  ac  (.round. 


Figure  11-6  -  Potential  Distribution  in  Interaction  Chamber  With  Target  at  -10  Volts 
and  all  Other  Electrodes  at  Ground. 


The  above  results  indicate  that  if  potential  differences  are  to  be  applied 
between  the  target  and  the  hemisphere,  the  following  precautions  should  be 
observed  in  order  to  minimize  the  distortion  of  the  angular  distribution  of 
particles  that  will  be  detected  at  the  hemisphere,  (a)  the  hemisphere  and 
the  outer  shell  of  the  chamber  should  be  held  at  the  same  potential  if  at  all 
possible,  (b)  a  fairly  good  approximation  to  a  radial  field  can  be  obtained 
by  minimizing  the  size  of  the  target. 

Since  with  an  applied  voltage  the  field  lines  that  run  out  from  the  target  to  the 
hemisphere  are  somewhat  distorted,  it  is  only  at  zero  applied  potential  to  the 
target  that  a  really  symmetric  case  is  obtained  This  indicates  that  the  best 
data  on  the  system  should  be  that  obtained  at  zero  accelerating  potential.  It 
is  difficult  to  see  how  the  results  obtained  in  Figure  11-4  can  result  from  a 
distortion  of  the  field  lines  such  as  are  shown  in  Figure  11-5  and  Figure  11-6. 
However,  it  does  indicate  that  possibly  the  best  results  are  those  in  the  top 
photograph  of  Figure  11-4.  The  angular  distributions  were  accordingly 
measured  with  no  potential  applied  to  the  target. 

C.  ANGULAR  DISTRIBUTION  MEASUREMENTS 

Working  at  low  laser  powers  so  that  the  data  were  free  of  noise,  we  have 

measured  the  angular  distribution  of  ions  emitted  from  the  target  in  the 

interaction  chamber.  Using  the  different  collectors  on  the  circumference 

of  the  hemisphere  and  taking  the  ratio  of  the  signal  from  a  detector  at  a 

particular  angle  to  the  signal  from  the  detector  along  the  normal  to  the  target 

surface  gives  the  angular  distribution  nor  malized  to  the  value  perpendicular 

to  the  surface  as  unity.  From  shot  to  shoi  at  constant  laser  power-  the  signal 

received  at  the  detector  in  the  direction  normal  to  the  surface  is  approxi- 

(5) 

mately  constant.  Preliminary  results  indicated  an  anisotropic  distribution 
with  more  material  emitted  at  75°  than  at  90  to  the  surface.  More  complete 
data  Iras  not  confirmed  this  preliminary  result.  The  m;  ximunr  amount  of 
emission  is  in  the  direction  nor  mal  to  the  surface  and  the  amount  falls  off  as 
the  angle  from  the  direction  to  the  normal  increases. 

W'e  have  obtained  the  angular  distribution  in  two  planes  perpendicular  to  each 
other  and  both  perpendicular  to  the  target  surface.  Since  the  laser  beam  slrrki 


the  target  at  an  oblique  angle,  the  direction  of  the  beam  and  the  normal  to  the 
target  define  a  plane,  which  is  one  of  the  planes  in  v  hich  we  made  measure¬ 
ments.  The  other  plane  is  perpendicular  to  the  first  and  also  contains  the 
normal  to  the  target  surface.  In  what  follows  we  shall  refer  to  these  two 
planes  respectively  as  the  plane  containing  the  beam  and  target  normal  and 
the  plane  perpendicular  to  thi  plane  containing  the  beam  and  target  normal 

The  results  are  shown  in  Figures  II- 7  and  11-8.  The  data  are  the  integrated 
emission  (i.  e.  area  under  the  curve)  for  a  particular  direction  divided  by 
the  ion  emission  from  the  collector  along  the  normal  to  the  target.  The 
normalization  reduces  scatter  due  to  shot  to  shot  variations  in  laser  power. 

Figure  11-7  shows  the  angular  distribution  of  ions  in  the  plane  perpendicular  to 

the  plane  containing  the  beam  and  the  target  normal.  The  angular  distribution 

2 

follows  approximately  a  cos  0  distribution  where  9  is  the  angle  between  the 
direction  of  observation  and  the  normal  to  the  surface,  and  is  symmetric  about 
0  =  0.  There  is  some  scatter  in  this  data  from  shot  to  shot,  but  the  general 
trend  of  the  data  is  well  apparent. 

Figure  11-8  shows  data  for  the  plane  containing  the  laser  beam  and  the  target 
normal.  The  amount  of  scatter  has  increased  considerably  and  the  angular- 
distribution  has  become  more  complicated.  We  regard  the  data  for  this 
plane  as  still  being  preliminary  and  subject  to  future  confirmation.  The  re¬ 
sults  in  Figure  II  -  8  indicate  an  asymmetry  about  the  normal  to  the  target 
surface.  The  laser  beam  enters  at  an  angle  equal  to  30°  from  the  target  nor¬ 
mal  through  the  hole  which  would  ordinarily  house  a  collector  in  that  direction 
so  that  no  data  were  obtained  at  that  angle.  At  30°  from  the  target  normal  on 
the  opposite  side  of  the  normal  in  the  same  plane,  the  laser  light  reflected 
from  the  target  struck  the  collector.  No  good  data  were  obtainable  from  this 
collector;  this  collector  gave  an  output  much  more  noisy  than  the  other 
collectors.  Apparently  the  reflected  laser  light  was  causing  particle  emission 
from  this  collector.  Accordingly,  no  data  are  given  for  this  angle  in  Figure  II -  8. 

From  Figure  Jl-8  it  appears  as  if  the  emission  with  a  component  along  the 
forward  direction  in  which  the  laser  beam  is  traveling  is  somewhat  less  than 
in  the  backward  direction,  at  least  at  the  angle  15°  from  the  normal  of  the  target 
surface.  At  angles  45°  and  greater  from  the  normal  to  the  target  in  both  the 
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forward  and  back  directions  the  emission  dropped  off  to  smaller  values,  so 
small  in  fact  at  60°  and  75°  that  it  was  virtually  unmeasurable  and  only  one 
po>nt  is  shown  in  Figure  11-8  for  these  angles.  It  appears  as  if  the  emission 
in  the  plane  is  concentrated  in  an  angular  spread  approximately  15°  from 
the  normal  to  the  surface  and  very  little  emer  es  at  larger  angles.  It  is 
unfortunate  that  the  data  at  30°  to  the  target  normal  is  not  obtainable  since 
this  would  fall  in  the  intermediate  region  between  relatively  large  and 
relatively  small  values. 

The  data  shown  in  Figures  11-7  and  11-8  were  both  taken  with  the  polari¬ 
zation  of  the  laser  beam  in  the  plane  containing  the  incident  laser  beam  and 
the  target  normal.  In  the  future  we  shall  take  data  with  the  laser  polarization 
perpendicular  to  this  plane.  It  is  not  immediately  apparent  lhat  this  change 
should  cause  a  difference  in  the  emission  but  it  is  a  point  worth  checking. 

Retarding  potential  measurements  have  also  been  carried  out.  Results 
indicate  that  the  signal  at  the  collectors  can  be  eliminated  by  application 
of  a  voltage  in  the  range  of  270  to  360  volts.  This  voltage  is  applied  on  the 
retarding  electrode  which  lies  between  the  screen  in  the  hemisphere  hole  and 
the  collecting  electrode.  The  screen  across  the  hemisphere  hole  shields 
ions  traveling  between  the  target  and  the  hemisphere  from  the  elfect  of  this 
potential  so  that  it  is  felt  only  by  ions  that  have  penetrated  through  the  screen. 
Figures  11-9  and  11-10  show  some  oscilloscope  traces  with  a  series  of  increasing 
applied  retarding  potentials.  The  results  are  es  .entially  in  agreement  with 
the  results  of  measurements  from  the  time  of  flight  spectrometer  for  the 
energy  of  ions  produced  under  similar  circumstances.  However,  we  note 
that  even  at  high  applied  voltages  on  the  retarding  grids  there  remains  some 
residual  signal  at  the  collector.  We  have  not  been  able  to  explain  this  signal. 

It  appears  to  be  independent  of  conditions  in  the  chamber.  At  270  volts  re 
tardmg  potential  the  smooth  shape  of  the  collector  lias  become  bumpy  and  the 
tail  of  the  pulse  has  apparently  been  cut  off. 
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Figure  II -9  -  Effect  of  Retarding  Potential.  Sweep  speed:  2|is/cm.  Top 
Traces  From  Collector  Along  Normal  to  Target,  0.  02V /cm. 
Bottom  Traces  From  Hemisphere,  lOV/cm. 

A.  0  volts  retarding  potential. 

B. .  90  volts  retarding  potential. 

C.  1  80  volts  retarding  potential. 
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Figure  11-10  -  Effect  of  Retarding  Potential.  Sweep  speed:  2|is/cm.  Top 
Traces  From  Collector  Along  Normal  to  Target,  0.  02V/ cm. 
Bottom  Traces  From  Hemisphere,  lOV/cm. 

D.  270  volts  retarding  potential. 

E.  360  volts  retarding  potential. 


The  angular  distribution  of  ion  energies  has  not  yet  been  investigated  by 
the  retarding  potential  technique.  We  shall  undertake  this  in  the  future. 

However,  some  information  is  available  from  the  same  data  that  were 
obtained  in  measurements  of  the  angular  distribution  of  particles.  A 
typical  oscilloscope  trace  of  the  output  from  two  different  collectors  is 
shown  in  Figure  II- 11.  In  A  the  upper  trace  is  the  emission  from  a 
collector  at  an  angle  of  15°  from  the  normal  to  the  target,  and  in  B  the 
upper  trace  is  from  a  collector  at  an  angle  of  75°  from  the  normal  and 
in  the  plane  containing  the  beam  and  the  normal.  Both  lower  traces  are 
from  the  collector  along  the  normal  to  the  target.  Both  the  amplitudes 
and  arrival  times  change  as  we  move  farther  away  from  the  normal.  The 
arrival  times  at  the  collectors  at  large  oblique  angles  are  later,  that  is, 
they  correspond  to  lower  ion  energies.  In  Figure  11-12  we  plot  the  relative 
energy,  as  derived  from  the  square  of  the  ratios  of  the  arrival  times  (i.  e. , 
the  time  at  which  the  peaks  of  the  pulses  occur)  as  a  function  of  angle  from 
the  normal  for  the  plane  perpendicular  to  the  plane  containing  the  beam  and 
the  normal.  In  Figure  11-13  we  show  the  same  results  for  the  plane  containing 
the  beam  and  the  target  normal.  These  preliminary  results  show  that  there 
is  a  variation  in  energy  of  the  ions  as  a  function  of  angle  of  emission. 

The  results  in  Figure  11-12  are  relatively  isotropic,  but  in  Figure  11-13 
the  energies  appear  to  be  maximum  in  a  cone  of  15°  around  the  normal  to 
the  target.  There  is  scatter  which  is  caused  by  variations  in  the  pulse 
shapes  from  one  shot  to  another. 

We  have  also  investigated  the  relative  numbers  of  particles  emitted  as  a  function 
of  laser  power.  Figure  11-14  shows  the  integrated  ion  current  from  the  collector 
normal  to  the  target  surface  as  a  function  of  laser  power  in  relative  units. 

There  is  a  considerable  amount  of  scatter  in  this  data  but  unmistakably  the 
general  trend  is  monotonically  increasing.  There  is  a  threshold  below  which 
no  emission  is  detectable. 
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Figure  II- 11  -  Ion  Arrival  Times.  Sweep  speeds:  5(is/cm. 

A.  Upper  Trace:  0.  OlV/cm  from  collector  at  15°  from  normal  to 
target  in  plane  containing  laser  beam  and  normal  to  target. 
Lower  trace:  0.  OlV/cm  from  collector  along  normal  to  target. 

B.  Upper  trace:  0.  005V/cm  from  collector  at  45°  from  normal  to 
target  in  plane  containing  laser  beam  and  normal  to  target. 
Lower  trace:  0.  02V/cm  from  collector  along  normal  to  target. 
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Figure  II- 13  -  Ion  Energy  Distribution  In  Plane  Containing  Laser  Beam 
And  Target  Normal 
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Figure  11-14  -  Ion  Emission  As  Function  Of  Laser  Power.  Integrated  Data  From 
Collector  Along  Normal  To  Target. 


HIGH  POWER  LASER 
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In  operating  the  new  large  ruby  laser  the  output  of  the  laser  was  raised  as 
high  as  200  megawatts.  However,  this  resulted  in  damage  to  the  ruby  rod. 
The  damage  appears  in  the  form  of  internal  voids  in  the  ruby  Some  of 
these  voids  which  scatter  light  in  a  transverse  direction  may  be  seen  in 
Figure  11-15,  Tne  photograph  was  taken  with  the  beam  of  a  helium- 
neon  laser  traveling  down  the  ruby  rod  axis.  A  possible  cause  of  this 
damage  was  filamentary  laser  action  in  the  rod  which  occurred  using  the 
Q- switching  dye  which  has  the  tendency  to  become  transparent  over  small 
areas  early  in  the  pulse.  This  leads  to  a  very  non-uniform  power  distri¬ 
bution  in  the  beam  and  local  flux  densities  much  higher  than  the  average 
value  can  occur.  This  means  that  in  small  filaments  in  the  rod  the  flux 
density  is  high  enough  to  cause  damage  even  though  the  average  flux  density 
over  the  face  of  the  rod  is  well  below  the  damage  threshold.  Photographs 
taken  in  a  Twyman- Green  interferometer  substantiate  this  view 

Figure  11-16  shows  a  Twyman-Green  interferogram  of  the  damaged  rod. 

The  rod  appears  as  the  roughly  elliptical  segment  near  the  center  of  the 
picture.  It  appears  this  way  because  of  the  Brewster  angle  end  faces. 
There  are  about  four  fringes  across  the  rod  diameter.  The  triangle  in 
'he  upper  left  is  a  background  reference.  The  damage  is  not  spread  out 
over  the  entire  surface,  but  is  localized  in  small  regions  across  the  cross- 
sectional  area  of  the  interferogram  substantiating  the  view  that  the  damage 
has  occurred  along  filaments  down  the  rod.  Moreover,  even  after  damage 
this  rod  was  operable  at  levels  of  approximately  50  megawatts.  This 
implied  that  it  is  only  a  portion  of  the  rod  which  underwent  the  laser  action 
at  flux  densities  in  excess  of  the  threshold  for  damage. 
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Figure  11-15  -  Scattering  Centers  In  Damaged  Ruby  Rod 
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Figure  11-16  -  Twyman-Green  Interferogram  Of  Damaged 
Ruby  Rod 
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A  new  six  Inch  long  ruby  rod  has  been  obtained  from  the  Linde  Corp.  In 
using  this  new  rod  we  shall  be  careful  to  monitor  the  laser  pulses  in  order 
to  avoid  exceeding  the  level  at  which  damage  to  the  rod  will  occur.  We 
have  set  up  new'  monitoring  procedures  in  which  we  have  calibrated  the  amount 
of  light  received  by  a  phototube  after  diffuse  reflection  from  an  MgO  block. 

By  keeping  the  entire  output  of  the  laser,  integrated  over  the  face,  below 
approximately  100  megaw'atts,  w'e  should  avoid  exceeding  the  level  at  which 
damage  occurs  in  the  filaments.  We  shall  also  use  a  single  cryptocyanine 
cell  instead  of  two  cells  within  the  optical  cavity  since  there  is  some  evidence 
that  filamenting  action  occurs  more  readily  when  there  is  a  cryptocyanine 
cell  at  each  end  of  the  laser  rod.  Because  of  the  noise  problems  in  the 
hemispherical  interaction  chamber  as  noted  above,  when  the  flux  density 
in  the  laser  beam  is  high,  we  shall  not  be  able  to  us  this  larger  laser 
immediately  to  extend  the  range  of  powers  in  our  measurements  as  wre  had 
originally  planned  to  do. 
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SECTION  III 


THEORETICAL  WORK 

A.  COMPARISON  OF  PULSE  SHAPES  WITH  FREE  EXPANSION 

We  have  compared  the  pulse  shapes  obtained  in  the  hemispherical  chamber 

to  the  results  predicted  if  the  pulses  arise  from  adiabatic  free  expansion 

of  a  high  temperature  plasma  produced  at  the  time  of  the  laser  pulse 

(5' 

immediately  in  front  of  the  target.  '  These  calculations  were  modified 

to  give  the  flux  as  a  function  of  time  at  a  position  9.  85  cm  from  the  target, 

a  distance  equal  to  the  radius  of  the  hemisphere.  Some  results  are  shown 

in  Figure  III- 1  for  the  asymptotic  case  of  adiabatic  free  expansion  of  a 

plasma  of  mass  23  ions  at  a  temperature  of  200,  000°  K.  This  is  compared 

to  the  observed  pulse  shape  obtained  from  the  collector  along  the  normal 

to  the  target  with  zero  voltage  applied  to  the  target.  From  earlier  measure- 
(1  2) 

ments  ’  we  deduced  on  the  time  of  flight  spectrometer  than  under  these 
circumstances  the  sodium  component  gave  the  dominant  contribution  to  the 
ion  emission  from  the  tungsten  target.  We  see  that  the  observed  pulse 
shapes  are  similar  in  shape  and  in  halfwidth.  The  predicted  value  from 
theory  has  a  somewhat  longer  tail,  but  it  may  be  that  the  tail  of  the 
experimentally  determined  pulse  is  lost  at  larger  time  values  because  of 
the  difficulty  in  determining  when  the  pulse  has  returned  to  the  base  line. 
Also,  the  earlier  rise  of  the  experimental  pulse  may  be  due  to  the  arrival 
of  a  light  component  such  as  hydrogen.  We  find  that  the  observed  pulse 
shape  is  approximately  half  way  between  the  calculated  values  for  the 
asymptotic  free  expansion  of  a  mass  23  plasma  at  100,  000°  K  and  200,  000°  K. 

These  results  are  in  agreement  with  the  model  in  which  the  blowoff  material 
(at  least  the  ions  in  it)  is  heated  to  — '20  ev  and  then  expands,  giving  the 
observed  pulse  shapes.  It  agrees  also  with  other  experimental  estimates 
of  the  temperature  of  the  blowoff  material,  based  on  optical  measurements, 

5 

which  appear  to  yield  values  of  the  order  of  10  degrees  under  a  variety 

(6  7  8) 

of  experimental  conditions.  ’  ' 
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Figure  lll-l  -  Pulse  Shapes  in  Interaction  Chamber 


B. 


PLASMA  FATING 


A  recurring  question  in  this  work  has  been  the  role  of  heating  the  blowoff 
material  by  incoming  laser  radiation.  At  relatively  high  laser  powers 
the  heating  rates  can  be  high  enough  and  the  particle  densities  large 
enough  that  a  considerable  amount  of  radiation  can  be  absorbed  in  the 
blowoff  material  and  the  production  of  high  energy  ions  can  be  expected. 

The  theoretical  treatment  appears  to  agree  reasonably  well  with  experi¬ 
mental  results  in  the  high  power  regime.  However,  ai  lower  powers,  around 
2 

50  megawatts/ cm  ,  there  has  been  a  problem  in  that  the  theoretical  calcu¬ 
lations  indicate  too  low  a  rate  of  energy  absorption  to  account  for  the  experi¬ 
mental  data.  We  analyzed  the  result  for  blowoff  material  consisting  of 
neutral  gases  such  as  we  observed  in  our  quadrupole  spectrometer,  and 
concluded  that  heating  by  absorption  of  the  laser  radiation  in  an  inverse 
Bremsstrahlung  process  in  the  neutral  gas  could  not  produce  the  observed 
results.  ^  Later  we  suggested  an  alternative  hypothesis,  namely  that 
the  absorption  of  radiation  takes  place  mainly  for  electrons  in  the  Coulomb 
field  of  alkali  metal  ions  such  as  we  have  observed  in  the  time  of  flight 
spectrometer.  The  electronic  energy  is  transferred  to  the  ions  via 
Coulomb  collisions  but  the  energy  is  coupled  relatively  inefficiently 
to  any  neutral  atoms  which  are  present.  This  model  results  in  an  increased 
rate  of  heating  at  lower  laser  powers  than  if  all  heavy  particles,  atoms  and 
ions,  shared  the  available  energy. 

Earlier  we  described  calculations  of  the  adiabatic  free  expansion  of  a 

(5) 

heated  gas  and  in  Section  II1-A  of  the  present  report  we  describe  a 
comparison  between  these  calculated  results  and  the  pulse  shape  observed 
in  our  interaction  chamber.  The  agreement  is  quite  good  for  expansion 
of  a  gas  heated  to  approximately  200,  000°  K.  This  leads  to  a  physically 
appealing  picture  of  heating  of  the  alkali  metal  ions  in  the  blowoff  material 
and  subsequent  free  expansion  of  the  alkali  ions.  However,  this  picture 
runs  into  difficulty  when  it  is  evaluated  quantitatively. 
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The  problem  has  been  approached  mathematically  from  two  points  of 
view.  First  there  is  a  computer  program  which  calculates  heating  and 
ionization  of  a  sodium  gas  heated  by  absorption  of  laser  radiation.  The 
second  approach  has  been  to  consider  the  heating  process  including 
equipartition  of  energy  in  the  ga:;  »nd  the  motion  and  spatial  extent  of 
the  material.  These  two  approaches  are  complementary.  Neither 
represents  the  full  problem,  but  each  gives  a  different  insight  into  it. 

Let  us  consider  the  computer  program  first.  In  an  ionized  gas  ai 
temperature  T  the  absorption  coefficient  for  radiation  of  frequency  v 
is  given  by  the  following  formula: 


3.  69  x  I08ne  [i  -  exp  (-hv/kT)]  E  Z2  N(Z)  log  A 

z 


(1) 


where  A  = 


2Z 


3e 

T 


k3T3  '  1/2 


rr  n 


ng  =  electron  density,  N(Z)  =  density  of  atoms  in  ionization  state  Z, 
and  the  sum  goes  over  all  ionization  states  present  in  the  gas. 


From  this  absorption  coefficient  we  get  the  rate  of  heating  of  the  gas, 
assuming  very  rapid  equipartition  of  the  absorbed  energy  between  the 
electrons  and  ions  and  neglecting  gas  expansion  during  the  heating. 


This  treatment  takes  into  account  the  fact  that  the  rate  of  heating  increases 
as  the  gas  becomes  more  ionized.  According  to  equation  (1),  as  one 
increases  the  level  of  ionization  from  one  to  two,  the  heating  rate  should 
increase  by  a  factor  of  four.  Therefore,  thermal  ionization  will  increase 
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the  heating  rate,  which  in  turn  gives  rise  to  more  ionization  in  a  process 

with  regenerative  feedback.  The  state  of  ionization  of  the  gas  was  con- 

(5) 

sidered  in  the  May,  1967  report  and  is  described  by  Equation  III- 4  of 
that  report.  This  represents  a  coupled  set  of  Saha  equations.  One  must 
first  calculate  the  partition  function  Z(T).  To  do  this  requires  knowledge 
of  the  excited  energy  levels  in  each  ionization  state  of  the  molecule.  These 
are  available  for  the  first  nine  ionization  states  of  sodium,  that  is,  for 

(9) 

sodium  ionized  eight  times  or  less,  in  the  NBS  tabulation  of  energy  levels. 

For  the  zeroth  ionization  state  of  sodium  there  are  over  170  excited  levels 
listed.  Since  each  excited  state  enters  with  a  weighting  factor  exp  (-E^/kT), 
where  T  is  the  temperature  and  E^  the  energy  of  the  nth  level,  this  means 
an  extremely  large  amount  of  computational  work  to  calculate  the  partition 
function  as  a  function  of  temperature,  and  thus  the  problem  is  suitable  only 
for  a  computer. 

A  Runge-Kutta  routine  was  written  to  calculate  the  heating  rate  of  sodium 
gas  absorbing  radiant  energy.  At  a  given  temperature,  with  the  populations 
of  the  various  states  known,  a  heating  rate  was  calculated  from  the 
absorption  given  by  equation  (1)  and  at  a  time  increment  later  the  new 
temperature  was  determined  according  to  the  Runge-Kutta  formalism. 

At  each  intermediate  step  in  the  Runge-Kutta  routine,  the  partition  functions 
for  the  intermediate  temperature  were  calculated  for  each  ionization  level 
of  sodium  in  a  subroutine  and  the  coupled  set  of  Saha  equations  was  solved 
to  adjust  the  populations  of  the  various  ionization  states  to  the  new  temp¬ 
erature.  The  electron  density  was  calculated  at  each  step,  and  the  change 

in  the  ionization  potential  due  to  the  Coulomb  interaction  as  described  in  the 
(5) 

May,  1967  report  was  also  calculated.  The  correction  to  the  ionization  po¬ 
tential  was  large  enough  so  that  for  each  ionization  state  of  sodium  enough 
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(9) 

energy  levels  were  given  in  the  NBS  tables  so  that  the  highest  value 
of  given  exceeded  the  corrected  ionization  potential.  Thus  we  did 
not  have  to  be  concerned  about  having  a  region  in  which  there  would  be 
a  large  number  of  unknown  energy  levels  slightly  below  the  ionization 
potential.  The  computer  routine  was  run  on  the  Honeywell  H-1800 
computer  for  various  choices  of  input  power  and  initial  sodium  density. 

The  starting  temperature  was  taken  to  be  5Q00°K,  roughly  characteristic, 
of  the  tungsten  surface  heated  by  the  laser,  but  results  are  not  critically 
dependent  on  the  initial  temperature.  At  5000°K  sodium  is  almost  entirely 
singly  ionized. 

We  had  expected  that  by  including  only  the  alkali  ions  in  the  system,  rapid 

and  effective  heating  would  occur,  particularly  when  we  allowed  for  the 

possibility  of  multiple  ionization  of  the  alkali  ions.  The  results  indicated 

that  heating  is  indeed  possible  at  lower  particle  densities  than  were 

indicated  by  the  results  of  the  similar  method  as  described  in  the 

(2) 

December..  1965  report.  Results  of  the  computer  programs  for  various 

initial  conditions  are  shown  in  Figure  III- 2.  The  laser  pulse  was  taken 

as  constant  in  time  with  duration  equal  to  30  nanoseconds.  The  various 

curves  represent  different  initial  conditions  of  particle  density  and  laser 

16  -  3 

flux  density.  Curve  a  represents  sodium  density  of  10  cm  and 

8  2  i7 

laser  flux  density  of  10  watts/cm  .  Curveb  is  for  10  cm  and 

8  2  is  -3  s  2 

10  watts/ cm  .  Curve  c  is  for  10  cm  and  10  watts/cm  .  Curve  d 

is  for  1018  cm  8  and  3  x  108  watts/ cm^.  Curve  e  is  for  10^ 8  cm  8 

9  2  “ 

and  10  watts/  cm  .  These  results  show  that  heating  can  occur  at  particle 

densities  at  least  an  order  of  magnitude  lower  than  indicated  by  the  earlier 

(2) 

model  in  which  neutral  gas  was  also  coupled  into  the  absorbing  system. 
However,  the  results  still  indicate  the  requirement  for  high  laser  powers 
and  high  particle  densities  to  cause  heating.  The  numbers  are  too  high  to  be 
realistic  for  our  system.  We  would  estimate,  from  measurements  of  the 
total  number  of  sodium  ions  in  the  time  of  flight  spectrometer  and  the  thermal 
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TEMPERATURE 

(°K) 
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velocity  of  the  ions  emitted  from  the  tungsten  surface  during  the  time  of 
the  laser  pulse,  that  a  density  around  10^  to  10^  sodium  ions/ cm'"*  would 
be  more  likely  in  our  experimental  conditions.  However,  ti  e  results  are 
interesting  in  that  they  indicate  the  conditions  under  which  plasma  heating 
can  occur.  It  is  worthwhile  to  consider  the  calculated  results  in  some 
detail. 

There  are  some  noteworthy  points  to  make  about  the  details  of  the  various 
computational  runs.  The  first  point  concerns  curve  e  in  Figure  III  -  2,  m 
which  conditions  for  a  runaway  heating  occur.  At  the  end  of  30  nano¬ 
seconds  the  computed  temperature  has  risen  above  400,  000°  K.  This  is 
probably  unrealistic  since,  as  the  material  gets  very  hot, it  would  expand 
rapidly  and  the  losses  would  limit  the  temperature  rise.  Indeed  it  is  very 
likely  that  the  approximately  constant  value  of  100,  000  to  200  000°  K  for 
the  temperature  of  the  blowoff  material,  as  observed  in  different  experi¬ 
mental  conditions,  may  indicate  the  level  at  which  losses  from  the  high  rate 
of  ovnan^ion  and  losses  by  radiation  become  large  enough  to  cut  off  the 
heating.  Thus,  the  potential  for  heating  by  this  mechanism  is  more  than 
adequate  to  account  for  the  experimentally  observed  energies  if  the  particle 
density  is  high. 

Another  point  worth  mentioning  about  curve  e  in  Figure  III -2  is  that  a 
high  enough  temperature  is  calculated  so  that  the  system  would  reach 
ionization  states  higher  than  eight  times  ionized  sodium.  These  higher 
ionization  states  are  cut  off  artificially  by  the  computer  routine  since  we 
do  not  have  the  available  input  data  for  the  energy  levels  of  the  higher 
ionized  states.  Thus  at  the  end  of  the  computation  we  are  in  a  physically 
unrealistic  situation  in  that  the  material  is  prevented  from  ionizing  to 
the  level  to  which  it  would  if  a  temperature  of  400,  000°  K  were  actually 
attained.  The  shapes  of  the  curve  are  also  interesting.  At  times  less 
than  about  eight  nanoseconds  the  curve  is  concave  downward.  This  occurs 
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because  of  the  form  of  equation  (1)  which  indicates  that  the  absorption 

1/2 

coefficient  decreases  with  increasing  temperature  because  of  the  T 
term  m  the  denominator,  so  long  as  the  ionization  state  does  not  change 
appreciably.  Wnen  the  temperature  reaches  about  30,  000"' K,  doubly 
lon'zed  sodium  becomes  an  appreciable  fraction  of  the  total.  When  this 
occurs,  the  heating  rate  increases  rapidly  and  the  curve  becomes  concave 
upward.  A  slight  upward  concavity  may  also  be  noted  just  at  the  end  of 
curve  d  where  a  temperature  of  approximately  35,  000°  K  and  a  considerable 
fraction  of  doubly  ionized  sodium  is  reached  near  the  end  of  the  laser 
puise. 

The  calculated  results  for  the  fractions  of  the  various  ionization  states 
are  of  interest.  Figu.  111.-3  shows  the  results  for  an  initial  density  of 
10  cm  and  a  flux  density  of  10  watts/cm*”,  i-  e. ,  for  the  same  case 
as  curve  c.  The  material  remains  almost  entirely  singly  ionized 
throughout  the  30  nanosecond  pulse  duration  '"•on  though  the  temperature 
is  calculated  to  rise  to  over  18,  000°  K.  The  material  is  practically  all 
singly  ionized,  but  the  fraction  of  un  ionized  material  rises  throughout 
the  time  duration  of  the  pulse.  This  is  somewhat  surprising  at  first 
in  that  we  would  expect  the  fraction  of  un- ionized  material  to  decline. 

The  increase  occurs  because  of  many  excited  states  of  un- ionized  sodium 
that  lie  at  relatively  low  levels,  below  35  000  wave  numbers.  Therefore, 
as  these  excited  states  become  accessible  as  the  temperature  rises  the 
partition  function  of  the  un  ionized  state  becomes  large  and  the  relative 
weight  of  this  state  causes  the  rise  as  shown  in  Figure  ID. -3.  Just  at  the 
end  of  the  laser  pulse  we  can  see  that  the  doubly  ionized  sodium  is  just 
coming  on  to  this  scale. 
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In  Figur  e  J!1  4  we  see  the  result  for  the  case  of  10  cm  and  a  flux 
8  2 

densttv  of  3  x  10  watts/cm“  i.  e.  for  the  same  case  as  curve  d.  At 
approx  matelv  25  nanoseconds  where  the  temperature  of  the  material 
has  risen  1o  approximately  30  000'  K,  doubly  ionized  sodium  has  reached 
approximately  ien  percent  of  the  total.  Triply  ionized  sodium  comes  on 
to  the  scale  at  a  low  level  near  the  end  of  fhe  pulse.  Un- ionized  sodium 
rises  to  a  broad  pcaK  because  of  tne  increase  ,jn  the  partition  function  as 
mentioned  before  ana  men  begms  to  decline  as  the  doubly  ionized  sodium 
begins  to  account  for  an  appreciable  fraction  of  the  material. 

In  Figure  111  5  w^  see  the  results  for  the  caie  in  which  we  have  runaway 

heating  to  a  very  high  level  i.  e.  for  the  case  of  an  initial  density  of 
18  3  m  2 

10  cm  and  lO"  w/cm  ,  the  same  as  curve  e.  This  curve  ts  cut  off  at 
approximately  23  nanoseconds  because  at  this  point  almost  all  the  sodium 
is  eight  times  ionized.  In  phvsicai  reality  we  would  at  this  point  have  nine 
times  ionized  sodium  commg  into  dominance  as  a  large  fraction  of  the  total 
number.  However  our  computer  program  has  no  input  data  to  account  for 
nine  times  ionized  sodium  and  tnercioie  cannot  realistically  account  for 
this  phenomenon.  Therefore  ve  do  not  present  the  results  at  times 
greater  than  23  nanoseconds.  We  see  that  the  behavior  is  qualitatively  as 
follows.  Singly  ionized  sodium  dominates  for  about  the  first  eight  nano¬ 
seconds.  At  11  nanoseconds  when  the  temperature  would  be  approximately' 

50  000°  K,  doubly  ionized  sodium  reaches  its  peak  and  begins  to  drop  off. 
Triply  ionized  sodium  reaches  a  peak  around  14  nanoseconds  when  the 
temperature  would  be  approximately  80  000°  K.  Four  times  ionized  sodium 
dominates  at  around  100  000"  K.  Near  150  000°K  (18  nanoseconds)  five 
times  ionized  sodium  is  largest.  At  200  000c  K  (20  nanoseeonds)  seven  times 
ionized  sodium  is  largest.  At  273,  000°  K  (  approximately  22  1/2  nanoseconds), 
eight  times  ionized  sodium  becomes  tne  largest  factor.  The  curves  re¬ 
present  the  beha1  tor  of  the  soluf’ons  of  the  coupled  set  of  Saha  equations 
111  4  of  the  May. 1067  report*^  as  the  tempei  ature  increases. 
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Figure  III- 4  -  Densities  for  Various  Ionized  States  of  Sodium  as  Calculated 

on  Plasma  Heating  Model  as  Functions  of  Time.  Initial  Sodium 
Density:  18  _3  8  2 

tO  cm  ;  laser  flux  density:  3  x  10  w/cm 
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The  results  have  the  following  implications.  If  we  do  indeed  heat  the 
bloivoff  material  to  approximately  100.  000°  K,  we  should  have  an  appreciable 
fraction  of  three,  four,  and  five  times  ionized  sodium.  We  have  not  detected 
these  species  in  our  experimental  system.  This  means  that  they  must  either 
recombine  very  fast  or  the  temperature  does  not  become  this  high.  It  is 
plausible  that  even  if  these  species  are  produced  they  could  recombine  so 
fast  that  they  would  not  be  observed  in  the  time  of  flight  spectrometer 
measurements.  They  should,  however  be  detectable  by  optical  emission 
lines.  Fawcett  for  example,  has  in  his  laser -produced  plasmas  seen 
lines  of  13  times  ionized  calcium  and  titanium. 

Thus,  on  the  basis  of  this  analysis,  we  conclude  that  the  conditions  for 

runaway  heating  to  temperatures  of  the  order  of  100,  000°  K  are  that  the 

18  3 

sodium  ion  density  should  be  approximately  10  cm  and  laser  flux 

9  2 

density  approximately  10  watts/ cm  .  We  could  postulate  the  possibilities 
of  nonuniformity  in  either  the  laser  power  or  in  the  sodium  density  as, 
for  example,  would  be  expected  in  the  tungsten.  This  would  give  rise  to 
a  local  very  high  concentration  of  sodium  ions.  However,  these  hypotheses 
are  rather  unsatisfactory  and  it  must  be  concluded  at  the  present  time  that 
there  are  grave  difficulties  in  the  quantitative  evaluation  of  this  model  in 
which  we  evaluate  the  possible  heating  of  the  blowoff  material. 

The  study  of  the  plasma  heating  problem  so  far  indicates  the  importance 
of  the  plasma  density.  Unfortunately  direct  information  on  the  plasma 
density  is  not  available  and  we  are  therefore  forced  to  leave  plasma  density 
as  an  independent  variable  or  construct  a  model  of  the  emission  process 
and  estimate  it.  In  what  follows  we  will  attempt  to  construct  a  useful 
model. 

In  the  treatment  described  before  we  neglected  many  important  points 
about  the  processes  occurring  in  the  blowoff  material.  In  particular  we 
neglected  any  motion  of  the  material  and  also  considered  the  equipartition 
of  energy  between  electrons  and  sodium  ions  to  be  instantaneous.  Including 
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these  factors  would  lower  the  calculated  heating  rate,  so  that  the  earlier 
results  give  an  overestimate.  We  will  attempt  to  introduce  the  particle 
densit}  m  a  more  natural  way.  instead  of  simply  allowing  it  to  be  an 
independent  variable. 

The  following  development  will  consider  some  of  these  points.  In  so 
doing  we  make  many  assumptions  which  will  appear  later.  Thus,  the 
treatment  is  limited  by  these  restrictions.  We  do  not  allow  the  realistic 
representation  of  the  increasing  ionization  with  increasing  temperature  as  we 
did  before,  but  consider  only  one  ionic  species.  The  result  is  subject  to  the 
stated  assumptions,  but  does  give  a  different  perspective  into  the  problem. 

We  assume  that  sodium  is  emitted  from  the  target  surface  singly  ionized. 

Any  atomic  or  molecular  material  (including  atomic  sodium)  will  not  be 
considered  as  part  of  the  system  because  of  the  long  time  scale  for  transfer 
of  energy  to  these  neutrals. 

The  process  considered  is  the  heating  of  the  electrons  by  inverse 
Bremsstrahlung  and  the  subsequent  transfer  of  energy  to  the  ions  by 
collisions.  Only  one  spatial  dimension  and  one  ion  species  are  considered. 

In  this  calculation  we  approximate  the  time  variation  in  laser  power  and 
target  temperature,  but  the  spatial  variations  are  ignored.  The  laser  pulse 
is  assumed  to  have  the  simple  form. 

t  st 

m 

t  sts2t 
m  m 

t  £  2t  (2) 

m 

where  P„  and  P  are  the  instantaneous  and  peak  laser  flux  densities  and 
i  m  r 

t  is  the  time  at  which  laser  power  is  maximum.  When  we  speak  of  "laser 


P„  =  P  t/t 

l  mm 


-  P  (2  -  t/t  ) 
m  m 


=  0 
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power"  we  normally  mean  average  (over  some  area)  peak  (in  time)  power. 
The  appropriate  value  of  P^  for  the  calculation  may  be  somewhat  higher 
than  this  average  peak  power. 


As  an  approximation  to  the  real  case,  it  will  be  assumed  that  the  instantaneous 

target  temperature  is  proportional  to  the  instantaneous  laser  power  density. 
(12) 

in  earlier  work  we  estimated  that  for  tungsten  targets  the  average  peak 

temperature  in  °K  is  approximately  60  times  the  average  peak  laser  flux 

2 

density  in  megawatts/ cm  for  a  30  ns  duration  laser  pulse.  The  target 
temperature  is  then  given  by 


0 


T 


S3 


60  x  10 


«  5  x  10'9  P,  (3) 

l 

where  0™  is  the  target  temperature  in  electron  volts  and  P  is  the  laser 
i  2  ** 

flux  density  in  watts/ cm  ,  and  we  take  e0T  to  be  of  the  order  of  kT,  where 

T  is  the  temperature  in  °K. 

By  making  assumptions  concerning  the  rate  of  ion  emission  and  the  velocity 
of  the  particles  an  expression  for  plasma  density  as  a  function  of  time  can 
be  obtained. 


We  assumed  that  the  X-component  of  velocity  of  the  ions  equals  the  average 
velocity  of  the  half- Maxwellian  distribution  corresponding  to  the  target 
temperature  at  the  time  of  emission. 
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2kT  \1^2 


Mtt 


(4) 


where  M  is  the  ion  mass.  From  this  point  on,  the  analysis  will  be  restricted 
to  the  first  half  of  the  laser  pulse.  Any  increase  in  velocity  due  to  the  plasma 
heating  during  the  first  half  of  the  laser  pulse  will  be  neglected. 
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At  time  t  <  t  ,  the  average  velocity  of  ions  emitted  t  seconds  earlier  »s 
found  by  substituting  equations  (2)  and  (3)  into  equation  (4). 

1 1/ 2 


v  (t,  t) 
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6  x  10' 5k 
2n  M 


m  (t 
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(5) 


and  the  position  of  these  ions  is  given  by: 


x  (t,  t) 


Vx  T 


(6) 


To  find  the  extent,  xmax-  of  the  plasma  at  the  time  t,  equation  (6)  may 
be  differentiated  with  respect  to  T  and  the  result  set  equal  to  zero.  By 
substituting  equation  (5)  into  equation  (8)  and  carrying  out  the  differentiation 
it  is  found  that  x  is  maximum  when  t  =  2t/ 3.  or 
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(7) 


Unfortunately  the  history  of  the  ion  emission  during  the  laser  pulse  is 
not  known.  A  crude  approximation  will  be  used  here  and  then  the  effects 
of  the  approximation  will  be  considered  at  the  end  of  the  analysis. 


It  will  be  assumed  that  the  ions  are  emitted  continuously  from  time  t  =  0 
at  a  rate  proportional  to  the  target  temperature.  Because  of  the  assumed 
time  dependence  of  the  target  temperature: 
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for  t  <  t  and  where  r  and  r  are  instantaneous  and  the  peak  ion  emission 
m  1  m 

flux  densities  respectively.  In  order  to  relate  the  ?  esult  to  quantities 
which  can  be  inferred  from  experimental  observations  it  is  assumed  that 


t 

m  P  di  ^  Nt/2 
i  T 

where  is  the  total  number  of  ions/ cm 
the  laser  pulse.  It  is  easily  shown  that. 


(9) 

emitted  from  the  target  during 


N, 


It 
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(10) 


Now  the  average  ion  density  will  be. 


n 
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Fi  dt 


(11) 


Upon  making  the  appropriate  substitutions,  for  M  =  23  times  the  proton  mass. 
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In  the  absence  of  laser  heating  the  average  electron  and  ion  temperatures 


0  and  0  would  be 
e  i 


0=0 
i  e 


(13) 


and  under  the  stated  assumptions 


0  =0 
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2  (0T)  t /  3  t 
T  max  m 


(14) 


For  convenience  we  have  used  the  average  densities  and  temperatures. 
Local  densities  and  temperatures  will  of  course  be  different. 


The  absorption  of  laser  energy  is  assumed  to  occur  through  the  inverse 

Bremsstrahlung  process  and  is  taken  to  be  described  by  the  expression 
(13) 

given  by  Spitzer.  The  absorption  coefficient  is  given  as 
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where  g^  is  a  correction  factor  "normally"  near  unity  and  will  be  taken 
as  such  in  this  calculation.  (Whether  or  not  the  situation  under  consideration 
is  "normal"  has  not  been  examined  in  detail.  )  This  differs  from  the  absorption 
coefficient  used  in  the  computer  program  in  that  here  we  consider  one  ionic 
species  to  be  dominant.  Inserting  the  numbers  and  symbols  appropriate 
to  this  analysis,  the  absorption  coefficient  is; 
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The  rate  of  increase  of  electron  temperature  due  to  absorption  is  given  by 


P„  k 

l  v 


Ions  will  be  heated  by  collisions  with  energetic  electrons  according  to  the 
equation 
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tgq  is  the  equipartition  time  for  energy  transfer'1  7  and  At  and  A^  are 
respectively  the  atomic  mass  numbers  of  the  ions  and  electrons.  For  the 
system  under  consideration,  A .  =  23,  Ag  =  1/1837, 

t  cs  7.  2  x  109  9  3/2  /Z2n.lnA  (20) 

eq  e  i 

(13) 

where  In  A  is  the  factor  as  defined  by  Spitzer  ,  approximately  equal  to 
10  under  a  wide  variety  of  conditi.ons. 

The  various  terms  can  be  collected  to  obtain  one  equation  for  electron  heating. 
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The  first  term  on  the  righf  is  due  to  the  rrnx’ng  of  the  hotter  ions  as  the 
target  heats  up  (obtained  by  d  fferenb  ating  equation  14),  the  second  term 
is  tne  laser  absorption,  the  third  represents  the  transfer  of  energy  to  the 
tons. 


The  ion  heating  may  be  expressed  by 
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Substitutions  from  rhe  preceding  analysis  can  be  made  into  the  above 
equations  with  the  result 
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These  two  equations  are  in  terms  of  given  parameters  P  .  t  ,  N™,  and 

*  m  m  T 

the  unknowns  9^,  0  and  Z  Toe  appearance  of  the  average  ionic  charge  Z 
considerably  complicates  the  solution  of  these  equations  since  it  can  be 
evaluated  only  by  the  use  ol  the  Saha  equations  as  discussed  previously.  ^ 
V  hile  the  simultaneous  solution  of  the  equations  of  this  analysis  could  be 
carried  out,  let  us  first  Iook  at  an  approximaiion. 
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The  maximum  electron  hea*  ng  r.i’e  would  be  near  t  -  t  .  An  approximation 

to  equation  <23)  at  t  =  t  would  be.  assuming  no  significant  heating  up  to  that 
n  m 

point: 
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This  result  is  interesting  in  that  the  laser  po*'  er  does  not  appear  explicitly. 
One  would  expect,  however,  that  the  value  cf  observed  in  an  experiment 
would  be  nigner  for  highe>"  laser  power  and  thus  higher  heating  rates  could 
be  expected. 

To  explain  the  observed  plasma  healing  on  the  basis  of  this  model,  d0  /dt 
would  have  to  have  a  value  of  the  order  of  10^  (1  eV/ nanosecond).  For  the 

-  Q 

experiment  t^  «  15  x  10  '.  )f  Z  «  1,  then  from  equation  (25), 
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N,j.  =s  2  x  10  sodium  ions/ cm  ^5) 


i.  e.  ,  this  is  the  value  of  N^.  for  which  heating  would  occur.  The  experi¬ 
mental  estimates  for  the  number  of  10ns  involved  in  the  heating  process 

9(2)  -  2  2 

are  about  10  10ns  If  they  are  emitted  from  an  area  as  small  as  10  cm 

11  2 

we  obtain  a  value  of  N.j,  about  10  lons/cm"  -  five  orders  of  magnitude  too 
low  to  explain  the  heating 

Of  the  various  approximations  used  in  the  analysis,  the  one  regarding  the 
ion  emission  rate  is  probably  the  least  satisfactory.  If  the  ion  emission 
were  more  like  a  step  function  or  a  delta  function  originating  at  a  certain 
target  temperature,  the  all  important  plasma  density  would  oe  higher.  A 
more  realistic  approximation  would  be  that  the  emission  is  an  exponential 
function  of  temperature.  Most  of  the  emission  would  then  occur  in  the  time 
interval  just  before  t  =  t  and  the  plasma  density  would  be  larger.  However 
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the  result  would  be  critically  dependent  on  an  activatton  energy  which  would 
not  be  well  known,  /t  is  unlikely  that  the  large  discrepancy  by  several 
orders  of  magnitude  could  be  resolved  in  this  fashion. 

In  order  to  compare  the  second  analysis  with  the  first  one  contained  m 

this  report  it  should  be  of  interest  to  calculate  what  plasma  densities 

are  predicted.  Taking  N„  =  2  x  10^.  P  =  1U^,  r  =  1.5  x  10 

then  equation  (12)  predicts  a  plasma  density  at  the  peak  of  the  laser  pulse 
IQ  3 

of  3  x  10  "  cm  .  This  is  not  ’nconsistent  wun  the  other  result.  Using  the 
experimental  value  of  =  10^  would  lead  to  =  1.5  x  10^^  cm 

The  equrpartition  time  is  also  of  interest.  For  the  conditions  assumed 

in  this  model,  and  for  the  experimentally  derived  values  as  in  the  preceding 

.  7 

paragraph  equation  (20)  predicts  an  equipartition  t’.me  of  approximately  10 
seconds.  The  observed  ion  heat  ng  would  require  an  equipartition  time  nearly 
two  orders  of  magnitude  shorter.  Again  the  plasma  density  is  the  dominant 
quantity,  a  nigher  densily  would  yield  a  shorter  equipartition  time. 

(5) 

Thus  on  the  basis  of  the  suggested  model  .  the  results  of  the  two  different 
approaches  although  the  assumptions  are  much  different,  are  qualitatively 
similar  m  suggesting  that  the  experimental  value  of  the  'on  density  is  too 
low  to  account  for  heating  of  the  ions  in  the  blowoff  material  to  temperatures 
of  the  order  of  20  eV  by  absorption  of  the  laser  radiation. 

C  PRESSURE  PULSES 

(5) 

In  the  May,  1967  report  we  described  results  of  calculations  using 

the  PUFF  computer  code  to  calculate  the  therinomechamcal  stress  profiles 

produced  in  tungsten  targets  irradiated  wnh  laser  radiation  in  30  nanosecond 

pulses.  The  result  typically  involved  propagation  of  short  sharp  stress 

8  2 

waves  with  ampi  tudes  of  the  order  of  6  x  10  dynes/ cm"  into  the  tungsten 
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target.  These  calculations  have  been  extended  to  nickel  and  carbon  target 
materials.  The  method  is  exactly  the  same  as  described  in  the  May,  1967 
report.  Some  results  are  shown  in  Figures  III- 6  and  III  -  7  for  the  stress 
as  a  function  of  depth  in  a  50  micron  thick  nickel  target  at  various  times 
after  the  start  of  the  laser  pulse.  Results  are  qualitatively  very  similar 

to  those  obtained  from  tungsten  targets.  The  peak  stress  is  approximately 

8  2  2 
6  x  10  dynes/ cm  for  a  total  energy  deposition  of  0.  3  cal/ cm  .  A  com¬ 
pression  wave  is  produced  in  the  early  stages  of  the  energy  deposition. 

The  compression  wave  travels  to  the  right  and  is  reflected  off  the  free 
back  surface  of  the  target  as  a  rarefaction  wave  traveling  to  the  left. 

A  gradual  attenuation  of  the  amplitude  of  the  wave  may  easily  be  seen. 

On  each  reflection  from  the  surface  of  the  target  there  is  180°  change  in 
the  phase  of  the  wave.  The  peak  stres..  produced  can  be  high,  of  the  order 
of  600  atmospheres,  but  this  is  less  than  the  yield  strength  of  the  metallic 
specimen  so  that  we  expect  no  permanent  deformation  of  the  material  to 
result. 

Similar  results  for  carbon  targets  are  shown  in  Figures  III-8  and  III- 9. 

In  these  figures  we  see  the  gradual  buildup  of  the  amplitude  of  the  original 
compression  wave,  and  the  start  of  travel  of  the  peak  amplitude  toward 
the  right.  The  peak  amplitude  in  this  case  is  reached  after  approximately 
half  a  nanosecond  at  a  value  of  approximately  0.  9  x  10  dynes/ cm  . 

Because  of  the  different  mechanical  properties  of  carbon,  the  amount  of 
computer  time  required  to  run  through  a  cycle  in  the  PUFF  code  and  to 
cover  a  given  increment  of  problem  time  becomes  unduly  large  and  this 
calculation  was  cut  off  before  the  pressure  wave  reached  the  rear  surface 
of  the  carbon.  We  would,  however,  expect  the  reflection  properties  of 
this  wave  to  be  exactly  analogous  to  those  of  the  other  materials. 
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We  see,  therefore,  that  we  expect  a  narrow  stress  pulse  propagating 
into  the  carbon  target.  The  magnitude  of  the  stress  pulse  is  somewhat 
smaller  for  carbon  than  for  the  metallic  elements  which  were  considered. 
Qualitatively,  for  materials  of  different  mechanical  properties  the  behavior 
of  the  stress  pulses  as  predicted  by  the  PUFF  code  is  very  Similar  but  the 
magnitude  can  differ  considerably. 

This  essentially  completes  our  consideration  of  the  pressure  pulses 
produced  in  various  targets  by  absorption  of  laser  radiation.  We  have 
no  plans  to  verify  these  computed  results  experimentally.  In  order  to 
carry  out  experimental  measurements,  the  procedure  would  be  to  attach 
pressure  transducers  to  the  rear  surface  of  the  specimen  to  be  irradiated. 
However,  since  the  frequencies  predicted  for  the  wave  in  this  analysis  are 
much  higher  than  the  resonant  frequencies  of  available  transducers,  we 
would  expect  the  experimental  problems  to  be  difficult  and  the  results  to 
require  a  considerable  amount  of  interpretation.  Therefore,  we  do  not  plan 
to  undertake  an  experimental  study  of  the  pressure  pulses  produced  by 
absorption  of  the  laser  radiation. 

D.  STEADY  STATE  HEATING 

Some  preliminary  calculations  have  been  made  of  the  temperature  that 

can  be  achieved  by  irradiation  of  different  surfaces  using  a  carbon  dioxide 

laser.  This  type  of  laser  is  very  desirable  in  laser- surface  interactions 

as  it  has  a  continuous  output  that  can  be  focused  to  produce  large  flux 

2 

densities  (of  the  order  of  about  0.  1  to  1.  0  megawatt/ cm  ).  This  type  of 
laser  can  also  be  Q- switched  repetitively  at  a  high  repetition  rate  and 
with  very  constant  pulse  shape.  Thus,  if  one  could  expect  to  achieve  heating 
of  the  plasma  from  the  surface  by  Q- switching  the  laser,  very  accurate 
measurements  of  plasma  temperature  as  well  as  particle  energies  and 
angular  distribution  could  be  made  using  this  laser. 
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In  the  following  calculations  it  is  assumed  that  we  have  a  thermally  thin 

disc  of  material  with  ?.  very  large  radius  that  is  heated  at  the  center  by 

the  focused  beam  from  a  carbon  dioxide  laser.  It  is  assumed,  further, 

that  the  flux  density  incident  on  the  disc  is  constant  out  to  a  radius 

where  ai^  is  the  gaussian  radius  of  the  focused  laser  beam.  This  results 

in  a  situation  in  which  we  effectively  have  a  line  source  of  radius  through 

the  disc  and  power  flows  out  of  this  source  into  the  disc.  We  have 

assumed  that  the  cooling  of  the  disc  takes  place  by  linear  convection 

from  both  sides.  We  have  neglected  radiation  because  of  the  difficulties 

in  treating  the  problem  with  non-linear  heat  transfer  terms.  Also 

neglected  are  the  phase  changes  in  the  material  which  is  considered 

as  an  idealized  solid  with  thermal  properties  independent  of  temperature. 

(14) 

The  problem  defined  above  is  treated  by  Schneider  The  temperature 

distribution  m  given  by. 
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Modified  Bessel  functions  of  the  second  kind 
Ambient  tmperature  in  °C 
Laser  power  output  in  watts 

Thermal  conductivity  of  material  in  watts/ cm°C 

Thickness  of  disc 

Radial  distance  from  center  in  cm 

Gaussian  radius  in  cm 

2 

Unit  surface  conductance  in  watts/ cm  0  K 
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The  equation  above  has  been  used  to  calculate  temperature  profiles  in 

different  materials  that  might  be  useful  in  laser-surface  interaction 

4  2 

studies,  assuming  a  unit  surface  conductance  of  6  x  10  watts/cm  °  K. 

The  materials  as  well  as  their  thermal  conductivities  are  listed  in 

Table  I.  The  temperature  profiles  obtained  for  20  mil.  thick  material, 

-  2 

assuming  a  gaussian  beam  radius  of  10  cm  and  10  watts  absorbed  in  the 
solid  are  shown  in  Figure  111-10.  The  temperature  profiles  shown  are 
those  corresponding  to  steady  state  conditions  i.  e. ,  the  temperature 
profile  that  will  be  obtained  after  irradiation  for  a  long  period  of  time. 

We  have  not  calculated  any  transient  temperature  distributions  yet. 

The  most  striking  feature  of  Figure  111-10  is  the  fact  that  the  maximum 
temperature  obtained  in  the  material  is  a  very  strong  function  of  the 
thermal  conductivity  of  the  material.  It  is  clear  from  the  figure  that 
for  a  10  watt  power  level  there  are  very  few  materials  from  which  one 
can  expect  considerable  evaporation  as  a  result  of  irradiation  with  a  laser. 
Furthermore,  the  figure  of  10  watts  absorbed  in  the  solid  is  much  more 
restrictive  than  it  would  appear  at  first  sight,  since  metals  have  extremely 
high  reflectivities  at  10.  6  microns.  Thus  in  order  to  obtain  10  watts 
absorbed  into  one  of  these  metals  a  laser  power  output  of  from  200  to 
500  watts  or  higher  may  be  required.  We  can  see  however  that  materials 
of  low  thermal  conductivity  such  as  graphite  and  plastics,  or  dielectrics 
in  general  can  be  raised  to  very  high  temperatures  with  a  modest  amount 
of  power.  In  this  case  10  watts  absorbed  in  the  material  means  slightly 
more  than  10  watts  output  from  the  laser  since  absorption  coefficients 
of  dielectrics  at  10.  6  microns  are,  in  general,  very  large  so  that  the  CO^ 
laser  provides  a  verv  efficient  heating  source  for  this  type  of  mater  ,il. 
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TABLE  1 


THERMAL  CONDUCTIVITY  OF  SELECTED  SOLIDS 


Material^ 

Glass 

Graphite 

Platinum 

Pyrolytic  Graphite  ( 1) 

Sodium 

Teflon 

Titanium 

Tungsten 


Thermal  Conductivity 
(watts  /  cm°K) 

0.  01 
1. 

0.  7 

0.  0093 
1.  35 
0. 0025 
0.  2 
1.  6 
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Figure  III- 10  -  Steady  State  Temperature  Profiles  Produced  in  Different 
Materials  Upon  Absorption  of  10  Watts  from  a  Gaussian 
Laser  Beam  Focused  to  a  Radius  of  10'  cm.  Calculations 
Assume  Power  Loss  to  Ambient  by  Linear  Convection. 
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SECY  *)N  IV 


SUMMARY 

The  work  described  in  this  report  deals  with  a  number  of  related  topics 
in  the  continuing  investigation  of  the  interaction  of  high  power  laser 
radiation  with  absorbing  surfaces.  This  report  treats  certain  selected 
aspects  of  this  study.  This  work  is  to  be  understood  in  context  with 
the  material  presented  in  earlier  reports  in  this  series  in  which  we 
described  measurements  on  the  particle  emission  produced  in  a  time  of 
flight  spectrometer  and  a  quadrupole  mass  spectrometer  (1,  2,  3,  4,  5). 
When  this  work  is  considered  together,  it  helps  lead  to  an  understanding 
of  the  nature  of  the  processes  occuring  when  the  high  power  laser  beam 
strikes  a  metallic  surface.  The  topics  described  here  are  extensions 
of  measurements  or  calculations  which  have  been  carried  out  previously 
as  part  of  this  project. 

The  specific  topics  treated  in  this  report  include  measurements  on  the 

interaction  chamber  in  which  a  number  of  collectors  are  located  around 

the  periphery  of  a  hemisphere  with  the  laser  illuminated  target  at  the 

center  of  the  hemisphere,  and  a  quantitative  theoretical  evaluation  of 

the  particle  emission  model  in  laser-surface  interactions  presented 

(5) 

earlier  in  which  the  high  energy  ions  are  rapidly  heated  by  absorp¬ 
tion  of  laser  light  in  the  inverse  Bremsstrahlung  process  with  very 
little  coupling  of  energy  to  the  neutral  molecules  in  the  blowoff  material. 

The  most  significant  results  in  this  report  period  are  the  following: 

(1)  Measurements  of  the  angular  distribution  of  the  ion  emission. 
The  ion  emission  is  maximum  in  the  direction  normal  to  the  target 
surface.  Tn  the  conditions  under  which  the  ions  were  measured, 
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they  consisted  mainly  of  alkali  ions.  In  the  plane  perpendicular 

to  the  plane  containing  the  laser  beam  and  the  target  normal,  the 

2 

angular  distribution  follows  approximately  a  cos  0  distribution 
where  0  is  the  angle  from  the  normal  to  the  target.  In  the  plane 
containing  the  laser  beam  and  the  target  normal,  preliminary 
results  indicate  a  more  complicated  distribution  with  most  of  the 
emission  being  contained  within  about  15°  of  the  normal  to  the 
target  surface. 

(2)  Measurements  of  the  angular  distribution  of  energies 

of  ions  emitted.  The  energies  of  ions  emitted  were  determined 
by  the  times  of  flight  to  the  hemisphere.  The  ion  energy  is 
maximum  normal  to  the  target  surface.  In  the  plane  perpendicular 
to  the  plane  containing  the  laser  beam  and  the  normal  to  the  target 
the  energies  are  nearly  isotropic.  In  the  plane  containing  the 
laser  beam  and  the  normal  to  the  target  surface,  the  energies 
fall  off  to  small  values  at  large  angles  to  the  normal.  The  absolute 
values  of  the  energies,  as  obtained  from  retarding  potential 
measurements,  indicate  agreement  with  earlier  results  obtained 
under  similar  experimental  conditions  on  the  time  of  flight 
spectrometer. 

(3)  Evidence  for  complicated  processes  in  the  interaction 
chamber  which  are  critically  dependent  on  the  voltages  applied 
to  the  system.  To  date  these  results  are  not  well  understood. 

(4)  An  indication  that  the  pulse  shapes  obtained  for  ion  emission 
in  the  interaction  chamber  agree  well  with  the  results  calculated 

on  a  model  in  which  the  blowoff  material  is  heated  to  approximately 
20  electron  volts  near  the  target  surface  and  then  is  allowed  to 
expand  in  an  adiabatic  free  expansion. 
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(5)  Quantitative  evaluation  of  the  model  in  which  the  heating  is 
produced  in  absorption  of  laser  radiation  by  electrons 

in  the  blowoff  material  in  the  inverse  Bremsstrahlung  process 

with  relatively  little  coupling  to  neutral  atoms  contained  in  the 

blowoff  material.  This  quantitative  evaluation  indicates  that 
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heating  occurs  at  particle  densities  around  10  cm  and  powers 
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of  the  order  of  10  watts/ cm  .  These  values  are  too  high  for 
our  experimental  condisions  and  it  appears  as  if  this  model  is 
not  adequate  to  explain  the  observed  phenomena. 

(6)  Calculation  of  thermomechanical  stress  profiles  in  the 
target  material  due  to  absorption  of  the  laser  light.  For  our 
experimental  conditions  the  main  feature  is  a  sharp  short  pressure 
spike  propagating  into  the  target.  The  shape  of  the  pressure  pulse 
is  relatively  independent  of  the  target  material,  but  the  maximum 
amplitude  is  dependent  on  the  material.  For  our  experimental 
conditions  the  maximum  amplitudes  do  not  exceed  the  yield  strength 
of  the  target  materials. 

(7)  Evaluation  of  use  of  a  continuous  CC^  laser  in  producing 
laser  effects.  From  temperature  profiles  calculated  on  a  steady 
state  basis  it  appears  as  if  the  COg  laser  may  be  a  valuable  tool  to 
use  in  an  extension  of  these  measurements,  particularly  for  target 
materials  where  the  free  electron  density  is  low.  For  such  materials 
the  COg  laser  can  provide  a  very  efficient  heating  source  with 
repetitive  pulses.  This  would  allow  measurement  of  the  particle 
emission  using  synchronous  detection  methods  and  possibly  elim-  • 
inate  some  of  the  noise  problems  experienced  using  a  ruby  laser 
with  a  low  duty  cycle. 

Future  plans  call  for  completing  and  extending  the  work  described 
in  this  report.  Specifically,  we  shall  complete  the  investigation 
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begun  on  the  interaction  chamber  and  carry  out  retarding  potential 
measurements  of  the  ion  energies  as  a  function  of  angle.  We 
shall  carry  out  the  experiment  to  determine  the  role  of  absorption 
of  laser  light  in  heating  the  blowoff  material  by  illuminating  this 
material  with  a  second  higher  power  laser  immediately  in  front 
of  the  target.  Optical  measurements  to  identify  the  ionic  species 
present  in  the  blowoff  material  through  their  emission  spectrum 
will  be  carried  out,  particularly  with  a  view  to  determining  the 
time  histories  of  the  various  species  present  and  to  determining 
if  the  ionic  species  are  compatible  v>.lh  temperatures  in  the 
blowoff  material  of  the  order  of  20  eV.  Finally,  we  will  investigate 
in  more  detail  the  application  of  a  repetitively  pulsed  (X^  laser 
to  this  type  of  measurement. 
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